Recent work of M.D. Johnston et al. has produced sufficient conditions on thestructure of a chemical reaction network which guarantee that the correspondingdiscrete state space system exhibits an extinction event. The conditionsconsist of a series of systems of equalities and inequalities on the edges of amodified reaction network called a domination-expanded reaction network. Inthis paper, we present a computational implementation of these conditionswritten in Python and apply the program on examples drawn from the biochemicalliterature, including a model of polyamine metabolism in mammals and a model ofthe pentose phosphate pathway in Trypanosoma brucei. We also run the program on458 models from the European Bioinformatics Institute's BioModels Database andreport our results.
展开▼